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kJSTRnCT 


A  program  for  the  routine  location  of  T-phase  sources  in  the  Pacific 
has  been  initiated.  Data  for  this  program  is  supplied  principally  by  the 
Pacific  Missile  Range  hydrophone  network.  Hydrophone  records  are  forwarded 
weekly  for  reading  at  Honolulu.  In  correlating  arrivals  records  from  all 
stations  are  viewed  together.  Correlation  is  established  by  similarity 
in  shape  and  level  and  by  locations  determined  roughly  from  arrival  time 
differences.  Arrival  times  and  power  levels  are  read  for  processing  hy 
an  IBM  70U0  computer. 

The  solution  for  location  and  origin  time  is  the  least-squares  fit 
to  all  hydrophone  arrivals  which  are  weighted  according  to  their  distrJHu 
tion  in  azimuth  and  their  distance  from  the  T-phase  source.  The  iterative 
solution  proceeds  on  the  assumption  of  zero  parallax  and  flat  earth  in 
the  vicinity  of  the  source. 

Velocity  is  derived  from  shot  calibrations  and  from  averaging  across 
a  contour  map  of  local  velocity.  For  the  computer  program  this  is 
expressed  as  a  power  aeries  in  latitude  and  longitude  for  each  station 
and  for  each  quadrant  of  the  Pacific. 

A  T-phase  strength  is  computed  fro-  readings  of  peak  power  level. 

This  strength  la  computed  for  a  distance  of  thirty  degrees  from  the 
source.  The  lower  half  and  the  upper  quarter  of  all  hydrophone  readings 
are  arbitrarily  rejected  from  computation  of  strength  as  being  influenced 
by  topographic  shadowing  or  containing  mistakes. 

Events  which  can  be  correlated  at  four  or  more  hydrophones  occur 
at  an  average  rate  of  about  one  per  hour.  Tabulated  solutions  will  be 
published  at  regular  Intervals. 


Studies  of  T  phases  in  the  Pacific  during  recent  years  have  shown 
that  many  naturally  occurring  seismic  events  were  defected  by  SOFAR 
hydrophone  +  hat  were  not  located  by  the  existing  seismograph,  network. 
Johnson  et_  al.  (1963)  reported  that  about  one-fifth  of  the  T-phase 
arrivals  at  a  single  hydrophone  station  could  be  correlated  with  earth¬ 
quake  epicenters  determined  from  seismograph  data  available  to  the 
Coast  and  Geodetic  Survey.  This  observation  implies  that  tne  use  of 
the  T  phase  in  seismic  source  location  might  significantly  reduce  the 
threshold  for  seismicity  studies  of  the  Pacific  and  its  borders. 


In  addition,  studies  of  the  T-phase  mechanism  require  a  detcnnh-i 
tion  of  the  T-phase  source,  as  distinct  from  the  earthquake  focus.  In 
order  to  satisfy  the  requirements  of  T-phase  research,  and  in  hope  of 
providing  a  supplementary  source  of  data  for  seismicity  studies,  & 
program  for  routine  location  of  T-phase  sources  in  the  Facific  was 
developed. 


Background 

A  statistical  method  for  earthquake  location  was  first  described 
by  Geiger  (1910).  With  the  advent  of  the  high-speed  computer,  appli¬ 
cation  of  the  basic  least-squares  technique  is  no  longer  laborious  and 
a  number  of  computer  programs  have  been  written  (Bolt,  I960;  Flinn,  I960; 
Gunst  and  Engdahl,  1962;  Steinhart,  19oli), 

Workers  in  underwater  sound  have  tended  to  favor  geometric  solutions, 
perhaps  because  the  number  of  recording  stations  is  usually  small  and 
also  possibly  because  a  lower  order  of  effort  is  required  to  develope 
the  computer  program.  Also,  these  programs  have  generally  dealt  with 
underwater  explosions  for  which  the  arrival  times  are  much  more  precisely 
defined  than  those  of  T  phases. 

The  reduced  precision  of  arrival  time  determination  in  the  present 
application,  together  with  the  high  incidence  of  events,  requires  that 
as  many  arrivals  as  possible  be  read  and  that  all  weighting  procedures 
be  incorporated  in  the  computer  program.  These  requirements  are 
satisfied  by  the  method  of  least  squares  which  is  employed  in  the  present 
program. 


Data  Collection 


Hydrophones  are  currently  boing  recorded  at  five  locations  in  the 
Pacific;  Eniwetok,  Wake,  Midway,  Oahu,  aid  California.  A  sixth  recording 
station  in  the  Aleutians  is  planned.  At  each  of  the  Pacific  Missile 
Range  stations  at  Eniwetok,  Wake,  and  Midway,  data  are  recorded  from  four 
widely- spaced  hydrophones  (-■  100  km).  The  California  record  is  from  two 
well-spaced  hydrophones .  Data  from  two  hydrophones  are  recorded  at  uahn, 
but  only  one  is  used  in  computations. 


Records  from  which 


arrival  times  and  T- phase  pof 


xevexs 


are  read 


are  roll-chart  displays 
charts  are  multi -c.ianne 


of  sound  power  level  versus  time  (Fig.  1).  The 
1,  from  those  stations  having  multiple  hydrophone 


sites . 


As  the  T  phase  is  strong  in  low  frequencies,  the  signal  is  passed 
through  a  15-cps  low-pass  filter.  This  filtering  improves  the  s_gnal- 
to-noise  ratio  by  10-20  deciuels.  Most  ship  noise  and  explosion  signals 
are  also  eliminated.  Because  of  the  low  frequencies,  a  low  stylus 
response  speed  is  used  in  order  to  give  a  smoother  trace  and  a  more 
meaningful  reading  of  power  level. 

To  reduce  the  bulk  of  records  as  well  as  paper  costs,  the  slowest 
available  chart  speed  is  used  (0.25  ran/sec).  This  practice  does  not 
compromise  reading  accuracy  as  the  T-pha3e  peak  is  typically  broad  and 
its  position  cannot  usually  be  defined  more  precisely  than  about  +_  3 
seconds. 

The  arrival  time  of  a  signal  is  defined  as  the  time  of  peak  power. 
Each  trace  is  kept  near  the  bottom  of  its  paper  channel  so  that  a  large 
amplitude  T  phase  will  not  lose  its  peak  against  the  stylus  stop. 
.Seasonal  variations  cf  sea-noise  background  may  make  it  necessary  to 
adjust  the  attenuators  so  that  the  traces  remain  near  the  bottom  of 
the  paper  channel  but  not  off  it. 

The  system  is  calibrated  daily  for  power  level  by  injecting  a 
10-cps  signal  in  place  of  the  sea  noise  at  the  input  to  the  hydrophone 
amplifier.  This  calibration  signal  is  varied  by  iC~db  steps  throughout 
the  full  dynamic  range  of  the  system.  It  appears  on  the  record  as 
shown  in  Figure  1.  This  calibration  signal  is  used  to  measure  the  power 
level  of  T  phases. 

On  each  chart,  one  margin  stylus  records  15-second  marks  from  a 
tuning  fork  chronometer  while  a  second  margin  stylus  is  actuated  by  the 
iiiO-and  600-cps  tones  of  radio  time  signals  (WWV,  WWVK,  and  JJY).  The 
radio  time  is  essential  to  detect  errors  of  the  local  chronometer. 

Chart  rolls  are  changed  nearly  simultaneously  at  all  stations 
according  to  a  twice-a-week  schedule.  They  are  mailed  to  Honolulu  at 
least  weekly. 

Data  Processing 


The  first  task  of  analysis  is  to  decide  which  T-phase  arrivals  at 
the  various  stations  are  from  the  same  earthquake.  This  is  made 
difficult  by  the  fact  that  there  are  usually  several  T  waves  crossing  the 
ocear  in  different  directions  at  the  same  time.  To  correlate  the  arrivals, 
all  records  from  a. particular  set  are  put  onto  a  large  chart  table  (Fig.  2). 
The  records  are  synchronized  and  cranked  across  the  table  one  hour  at  a 
time. 


? 


Correlation  is  done  by:  (1)  sc -lining  the  synchronous  chart  display 
for  T  phases  of  similar  size  and  shape  within  a  possible  time  range;  (2) 
comparing  the  bearings  at  each  station  as  given  by  the  phone -to- phone 
differences  in  arrival  time :  and  f3)  comparing  the  station-to-station 
differences  of  arrival  time  with  the  Pacific  Ocean  chart  above  the  table 
which  shows  position  lines  for  given  time  difference.  Distances  from 
suspected  soiree  areas  are  checked  on  a  81.28  cm  globe  with  a  tape- 
calibrated  in  SO  FAR  travel  time. 


When  a  signal  is  found  which  correlates  at  four  or  more  hydrophones 
arid  at  two  or  more  stations,  it  is  assigned  a  two-digit-  serial  number. 
This  number  is  written  on  the  record  over  each  hydrophone  trace  which 
shows  the  signal  clearly.  The  serial  numbers  and  approximate  source 
areas  are  logged  for  future  reference. 

After  the  correlations  have  been  completed,  the  records  are  ready 
for  digitizing  on  an  automatic  char  reader.  This  chart  reader  is 
connected  to  a  card- punching  machine,  which  translates  the  data  to  IBM 
cards.  The  automatic  chart  reader  handles  one  roll  at  a  time  and  feeds 
the  following  information  to  the  card-punching  machine;  serial  number, 
power  level,  date,  arrival  time,  and  hydrophone  number.  It  is  only 
necessary  for  th„  operator  to  align  cross  hairs,  turn  switches,  and 
press  a  button.  One  card  is  punched  for  each  T  phase  at  each  hydrophone. 


One  of  the  "cross  hairs"  which  the  operator  aligns  when  measuring 
power  level  is  actually  a  pencilled  curve  on  a  plastic  sheet  which  he 


constructs  from  the  power-level  calibration.  Each  time  he  encounters  a 
calibration  in  the  chart  roll,  he  checks  his  curve  to  see  whether  it 
agrees  with  the  new  calibration.  If  it  does  not,  this  indicates  that  the 
system  characteristics  have  changed,  perhaps  due  to  recorder  sensitivity 


or  centering  adjustments,  or  to  attenuator 


As  this  is  precisely 


the  purpose  of  the  calibration,  the  operator  reconstructs  the  curve  to 


agree  with  the  new  characteristics. 


After  the  data  from  a  complete  set  of  chart  rolls  have  been  punched 
on  cards,  the  cards  are  machine  sorted  by  serial  number  and  by  arrival 
time  within  serial  number.  They  are  then  ready  for  the  electro  .ic 
digital  computer  which  solves  for  source  coordinates  and  origin  times. 


Computational  Method 


The  machire  program  first  checks  arrival  time  differences  and  the 
number  of  readings  in  eacv  set  to  detect  gross  errors  in  chart  reading. 
It  then  proceeds  to  select  four  hydrophones  for  use  in  a  geometric  fix. 
This  geometric  fix  serves  as  the  smarting  point  for  an  iterative  least- 
squares  estimate  of  the  source  which  uses  all  arrival  times. 


It  is  commonly  thought  that  three  detectors  are  sufficient  to 
obtain  a  fix  where  velocity  is  assumed.  However,  three  hydrophones 
give  two  solutions  (a  hyperbola  on  the  surface  of  a  sphere  is  a  closed 
curve)  and  a  fourth  hydrophone  is  still  required  to  resolve  the  ambiguity. 


The  first  hydrophone  selected  for  the  four-hydrophone  fit  is  the 
one  with  the  earliest  arrival  time.  The  second  is  tho  earliest  hydro¬ 
phone  among  the  remaining  recording  stations.  If  the  signal  was  recorded 
at  more  than  two  stations,  the  process  continues;  if  not,  the  third 
hydrophone  is  the  sscond  earliest  at  tho  station  recording  the  earliest 
signal.  If  only  one  hydrophone  was  recorded  at  that  station,  the  choice 
passes  to  the  succeeding  station,  etc.  The  objective  is  to  select  four 
hydrepnones  as  wide-spread  as  possible  but  as  near  the  source  is  possible 
where  alternatives  exist. 


k 


Fig.  3 


Although  it  is  possible  to  obtain 
latitude,  longitude,  origin  time,  and 
average  velocity  from  four  hydrophone 
arrivals  (Frosch  et  a]..,  1961),  the 
purpose  of  the  solttion  requires  only 
that  an  approximate  latitude  and 
longitude  be  obtained.  Therefore,  for 
this  preliminary  solution,  an  average 
SOFAR  velocity  of  2.32  x  10“u  geocentric 
radians  per  second  is  assumed  to  apply 
throughout  the  Pacific. 

Four-Station  Fix 


Consider  hydrophones  designated 
0,  1,  2,  and  3>  and  Increasing  distances 
6,  Ax,  As,  and  63,  respectively,  from 
the  source  at  k.  Distances  from  hydro¬ 
phone  0  to  the  other  hydrophones  are 
Bj_,  B2,  arid  B3,  respectively,  and  angles 
0,  »  ,  and  3  ,  are  as  shown  in  Figure  3. 

The  basic  equations  are  the  set 


H  •  T 


i 


for  origin  time  B,  arrival  time  T^, 
average  velocity  c,  and  i  corresponding 
-  ■)  hydrophones  0  through  3* 


Eliminating  H  between  equations  0  and  1  and  multiplying  by  c, 


cCT^  -  T)  ♦  6  •  «  eos’keos  cos 

Taking  the  cosine,  expanding  the  left-hand  sine. 


6  *  sin  B-i  sin  6  cos  0) 
dividing  by  sin  6  ,  and 


(1) 


-h- 


.ng  for  cot 


0»  S 


-Ves 


cot  A  = 


sin  Pi  +  si n  3-j  cos  0 
cos  D-|  -  cos  3^ 


where 


>x  -  c(Tx  -  T). 


By  parallel  construction, 


cot  A 


sin  D-  ♦  sin  B2  cos  (0  *  <*} 

cos  D,->  -  cos  B., 

^  2 


Eliminating  A  , 


£~  sin  D-,  +  E-s  sin  Bn  cos  0  ■  E,  sin  D0  +  EL  sin  3„  cos  (0  +  <*) 

t  1  i  t  i  ^ 


where 


*  cos  D-^  -  cos  B., 


Expanding  cos  {'6  +  <*)  and  collecting  terms, 


cos  0  ( E,-.  sin  B,,  -  E,  sin  3  cos  01)  +  sin  0(E  sin  B  sin  ») 

*  1-1  2  1  2 


E^  sin  D0  -  sin 


By  parallel  construction, 


cos  0  (E^  sin 


E  sin  B  cos 
1  3 


)  +  sin  0(E  sin  B  sin  3 )  * 
1  3 


Ei  sin  -  E-j  sin  D, 
1  >  *  l 


These  equations  may  be  solved  for  cos  0  and  sin  0  which  may  in  turn  oe 
combined  to  yield  a  single  equation  in  velocity.  As  the  purpose  is  to 
achieve  3  starting  point  for  the  iterative  least-squares  estimate,  the 
solution  for  velocity  does  not  seem  profitable.  Rather,  an  approximate 


[oiiaea  irom  the  "uot-iej 


sin  0  and  -?«s  0,  yields  good  resolution 


•hroughout  the  circle.  The  quadrant  of  0  is  determined  from  the 
algehraic  signs  of  the  numerator  and  denominator. 

Distance,  A  ,  may  be  obtained  from  (1)  or  its  two  parallel 
equa  ions,  each  of  which  involves  a  different  paii  of  hydrophones. 

In  cases  where  two  hydrophones  are  close  together  in  azimuth  from  the 
source,  the  small  velocity  error  may  cause  a  large  distance  error.  To 
avo  d  the  selection  of  such  a  pair  of  hydrophones,  distance  is  computed 
-  rcm  i  jj,  .he  three  equations  and  the  median  value  is  chosen. 

Veloc ity 


-^ihal  stxution,  velocity  is  dete  mined  from  a  function  of 
>  rop  one  i  and  p«xnt  k.  This  function  is  th>  least-squares  fit  to 
average  velocities  determined  from  such  ti  n.  aaliorations  as  are  avail¬ 
able  and  from  velocities,  averaged  across  a  contour  chart  of  local 
velocity.  Figure  L,  (Johnson  and  Norris,  196k),  to  selected  j  nts  in 
seismic alxy  active  areas  oi  tne  Pacific.  The  function  is  a  power 
series  oi  degree  2  in  latitude  L  and  longitude  ^  of  the  form 


for  h  pc in zs  to  which  velocity  has  been  measured,  p  *  0,  1,  2  and 

q  *  o,  i,  2. 


A  separate  set  of  coefficients  was.  obtained  for  each  quadrant  of 
the  Pacific  and  a  separate  set  of  functions  was  derived  for  each 
recording  station. 

Calculation  of  velocity  in  geocentric  radians  for  specific  recording 
stations  obviates  the  consideration  of  ellipticity  corrections. 


The  least-squa-es  estimate  of  the  T-phase  source  is  that  latitude 
L  a-id  longitude  \  wnich  yields  the  minimum  estimate  of  the  variance 
*u  origin  times  to  the  set  cf  equations 


This  definition  is  completely  equivalent  to  that  customarily  used  in 
determining  earthquake  epicenters— that  the  epicenter  is  given  by  that 
origin  time,  latitude,  and  longitude  for  which  the  sum  of  the  squares 
of  the  station  residuals  (observed  minus  computed  arrival  time)  is 
minimum.  The  present-  definition  has  the  advantage,  however,  that 
it  is  not  necessary  to  initially  assume  a  trial  origin  time  as  thx 
quantity  is  supplied  by  the  weighted  mean  of  the  computed  origan  times. 

Weighting  would  be  best  achieved  from  an  analysis  of  the  computed 
deviations  in  origin  time  were  the  number  of  recording  hydrophones 
statistically  significant  and  uniformly  spread  in  azimuth.  However,  in 
the  present  hydrophone  network,  a  T  phase  for  which  a  source  is  computed 
is  well  recorded  at  an  average  of  six  hyd“ophones.  The  number  of 
independent  recording  stations  is  ever  less.  Also,  it  is  the  nature 
of  T-phase  observation  that  hydrophones  are  confined  to  a  small  sector 
of  azimuth  with  respect  to  the  T-phase  source.  Therefore,  it  seems 
appropriate  for  the  present  problem  to  weight  according  to  consideration-- 
of  geometry. 

Hie  distance,  A  f  is  subject  to  uncertainties  c  5  T  and  A  6  c/c 
where  6  T  arid  6  c  are  uncertainties  in  arrival  time  and  velocity, 
respectively.  Combining  these  according  to  the  rule  for  addition  of 
errors,  the  uncertainty  in  A  is 


Si  =,/  (c  6  T)2  +  (A  —  f 

c 

The  uncertainty  of  arrival  time  varies  considerably  from  event 
to  event  according  to  the  sharpness  of  the  peak  of  power  level.  To 
assess  this  quantity  objectively  for  each  arrival  does  not  appear 
practical  considering  the  high  incidence  of  events  and  the.  available 
labor  force.  A  value  of  three  seconds  is  considered  representative  and 
is  applied  to  all  arrivals. 

Similarly,  the  uncertainty  in  velocity  varies  with  the  path  as 
certain  areas  are  better  calibrated  than  others  and  the  functions  from 
which  velocities  are  computed  give  varying  fits  to  the  calibrations. 
Here,  a  value  of  0.0005  is  considered  representative  of  6 c/c  and  is 
applied  to  all  paths. 

If  a  pair  of  hydrophones  is  considered,  the  uncertainty  xn  the 
hyperbola  of  position  of  the  source  depends  both  on  the  uncertainty  of 
distances  to  the  two  hydrophones  and  on  their  difference  in  azimuth. 

In  Figure  5,  with  hydrophones  i  arid  m,  the  hyperbola  of  position  is 
shown  passing  through  the  trial  source  location  k.  V  ^  is  the 


-9- 


differenc#  in  asiauth  fro*  k  of  the  hydrophone  p*ir.  The  uncertainties 
la  distance  fora  a  parallelogram  with  side  ^  .,/sin  v  and  /sin  4 

The  projection  of  this  parallelogram  upon  the  normal  to  loe  hyperboTa  is 
ft  Measure  of  the  uncertainty  in  the  position  of  the  hyperbola. 


Pig.  5 


If  the  hyperbola  is  approximated  by  the  bisector  of 


«ml  =iAkl 


_  CO;,  X&l  +  cos 

sin  Ymi  2  sin  Yni  2 


mi 


2U 


iLr 


ki 


6h, 


sm 


sin 


"^iai 


the  weighting  function  is  taken  aa 


-7 


ra-l 


2U. 


sin£ 


iii 


m-i 


for  arrival  times  read  at  n  hydrophones. 


-10- 


Leas  t-Squarea  Sstlmate 

In  obtaining  the  least-squares  estimate  of  source  location,  the 
formula  for  distance  Is  linearized  by  assuming  zero  parallax  and  flat 
earth  in  the  vicinity  of  trial  source  point  V .  As  shovn  in  figure  6, 


Fig.  6 


•  Akl  ♦  (  x  -  \k)  cos  1^  sin  -(L  -  1^)  cos  ©^  (3) 


for  azimuth  ©kl.  Steinhart  (196L),  in  solving  a  related  problem, 
arrived  at  an  equivalent  linearization  through  a  Taylor  expansion  of 
the  exact  formula.  The  coefficients  of  the  latitude  and  longitude 
differs nees  can  alto  be  obtained  by  differentiating  the  exact  formula 
(Hodgson,  1937). 


Substitution  of  (3)  into  (2)  leads  to  the  least-squares  estimates 
of  \0  and  I* 


\w^  3^  )  V,  FI  -  v±  FG  )  GI 
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The  new  coordinates  of  the  source,  L  and  X  ,  are  compared  with 
the  previous  estimate,  and  Xi  .  If  either  differs  by  0.1  degrees  o 

more,  then  L,_  and  X,  are  replaced  by  Lc  and  and  the  iteration 

continues,  if  both  differ  by  less  than  0,1  degrees,  the  position  is 
accepted,  distances  are  computed  and  an  estimate  of  origin  time 
is  obtained  as 


H 


LQ  is  converted  to  geodetic  latitude  before  print-out. 


the  machine  ti 
the  recording 
convergence  is 


ory  estimates  are  usually  obtained  in  three  iterations  and 
me  per  event  is  about  five  seconds.  However  in  cases  where 
stations  are  confined  to  a  small  sector  of  azimuth,  the 
not  adequate  for  the  flat-earth  approximation  and  the 


iterative  estimates  may  oscillate.  If  the  oscillations  are  of  decreasing 
amplitude  the  orogram  is  allowed  to  continue.  If  the  machine  program 
detects  increasing  amplitudes,  then  the  last  two  estimates  are  averaged 
and  iteration  proceeds  from  the  midpoint.  If  oscillation  of  inert  asing 
amplitude  again  takes  place,  then  the  average  of  the  last  two  estimates 
is  printed  with  lower  precision  indicated  by  one  le3s  significant  figure 
in  latitude  and  longitude.  The  averaging  process  is  also  used  to  accel- 
lerate  convergence  if  oscillations  with  decreasing  amplitude  require 
more  than  six  iterations. 


Measurement  of  Precision 


If  we  consider  the  variance  of  origin  times,  as  a  quantity  whioh 
varies  continuously  with  L  and  X,  then  L0,  XQ  is  that  point  at  which 
V  is  a  minimum  (  "  v  *  0).  I*  can  be  shewn  that  any  vertical  section 
through  the  variance  surface  is  parabolic  in  the  vicinity  of  L0,  X0, 

The  shape  of  such  a  section,  then  is  defined  in  the  vicinity  of  L0,  XQ, 
by  two  quantities,  V,  and  its  second  derivative.  The  variance  surface 
at  Lq,  Xq,  is  characterized  by  V  and  the  LaplAcisn  of  V 


il  .  sec2 
dL2 
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d2V 
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t:n 


(n-1) 


w« 


where  G  and  F  are  as  defined  in  (U)  and  (5). 

It  will  be  noticed  that  d is  independent  of  arrival  times.  It 
serves  as  a  measure  of  the  accuracy  of  the  geometry  of  the  fix^  the 
position  being  more  sharply  determined  for  larger  values  of  7<-V.  Bolt 
(i960)  recognized  the  usefulness  of  such  a  measure.  Examples  of  calcu¬ 
lations  giving  high  and  low  values  of  are:  (a)  An  earthquake 

T-phase  source  in  the  Marcus-Necker  Seamount  Chain  at  18.?  N,  176.8  E, 
Qlh  uiim  36s  GMT  on  2?  January  1?65  was  recorded  on  four  hydrophones  at 
Enivetok,  three  at  Wake,  and  three  at  Midway.  d^y  was  calculated  as 
3ii,3  seconds VmilliradianS  Independent  calculations  using  P  waves 
gave  an  epicenter  with  the  same  coordinates  and  one  second  earlier  in 
origin  time,  (b)  A  T-phase  source  was  computed  for  a  Solomon  Island 
earthquake  at  11.2  S,  159.1  E,  20h  OLtm  iris  GMT  on  22  August  196ii,  It 
was  recorded  on  one  hydrophone  at  Eniwetok,  two  at  Wake,  and  one  at 
Midway.  d^Y  was  1.0.  The  C  &  GS  epicenter  was  at  10.1  S,  161.7  E, 

Oh  7m  51. 7s.  Figure  7  shows  the  geometry  of  both  fixes. 


e  to  occur  when  d^V  is  less  than  2. 


Oscillation  is  . 


-Phase  Stre 


neasursments  oi  power  level  of  the  T-pftase  pe&x  at  the  various 
hydrophones  are  used  to  compute  a  T-phase  strength — a  concept  similar  to 
eart  quake  magnitude-  The  T-phase  strength  i_  intended  to  represent  the 
peaic  power  level,  in  decibels  relative  to  0.1  microbar,  of  the  T  phase 
at  a  distance  of  30°  from  the  source  on  an  unobstructed  path,  and  as  seen 
through  a  system  with  a  frequency  response  equivalent  to  that  of  the 
Pacific  Missile  Range  stations.  The  effect  01  ^he  PMR  system  response 
is  to  shift  the  peak  of  the  T-phase  power  spectrum  to  about  10  cps.  The 
0.1  mierobar  reference  was  chosen  so  that  all  power  level  measurements 
would  be  positive.  The  30°  distance  was  chosen  as  the  length  of  a 
typical  T-phase  rath  from  source  to  hydrophone.  Such  adjustments,  then, 
as  are  necessary  to  compute  the  power  level  at-  the  standard  jCp 
are  minimized. 

Reduction  of  each  reading  to  the  standard  distance  is  accomplished 
by  the  formula 


Si  *  Pi  +  10  log  2  sin  A  4  +  11-1  (  - ~) 

for  standard  power  level  S,  and  measured  powsr  level  P.  The  second 
term  on  the  right,  hand  side  of  the  equation  accounts  for  spatial 
spreading  over  the  earth's  surface  while  the  third  term  accounts  for 
losses  at  a  rate  of  1.6  db  per  megayard  as  reported  by  Urick  (1961i) 
for  a  frequency  of  10  cps.  For  the  computer  solution,  A  .  is  in  radians 
and  the  loss  coefficient  is  shown  in  appropriate  units  (db/ radian). 

A  point  source,  such  as  an  underwater  explosion,  would  require  an 
additional  term  to  account  for  spreading  of  the  signal  in  time  (Urick, 
ibid).  The  earthquake  T-phase  source,  however,  is  spread  in  both  time 
and  space,  and  the  T-phase  peak  is  typically  broader  than  the  rise 
time  of  the  explosion  S0FAR  signal.  Time  spreading,  then,  does  not 
affect,  the  measurement  of  peak  power  of  the  T  phase. 


Depending  on  the  location  of  a  particular  source,  arrivals  at 
certain  hydrophones,  although  readable,  are  along  paths  partially 
obstructed  by  topographical  features.  The  power  levels  at  those  hydro¬ 
phones  are  consequently  low.  Also,  the  incidence  of  mistakes  leading  to 
inordinately  high  readings  of  power  level  is  not  negligible.  As 
the  volume  of  data  makes  the  subjective  investigation  of  each  case 
prohibitive,  and  as  the  minimum  number  of  hydrophones  used  in  determining 
a  T-phase  source  is  four,  it  seems  appropriate,  in  the  case  of  the  four- 
station  fix,  to  arbitrarily  eliminate  the  lowest  two  and  the  highest 
standard  power  levsls  and  to  accept  the  next  highest  value  as  the  T-phase 
strength. 


To  extend  this  evaluation 
formula  has  been  devised 


fixes  from  a  hydro phones  the  following 
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where  brackets  signify  that  the  next  smaller  integer  value  is  to  be 
used  if  the  value  within  the  brackets  is  net  an  integer.  is  ordered 
in  ascending  values.  S  is  the  T- phase  strength. 

It  is  noted  that  the  evaluation  of  T-phase  strength  does  not 
consider  the  degree  of  coupling  of  body  waves  to  T  waves  at  the  source, 
azimuthal  radiation  patterns  at  the  source,  nor  the  relative  efficiencie 
of  the  wave  guides  to  the  various  hydrophones. 


Simplified 


-  J.UW 


u  ai-i 


Read  constants 


differences 

>0.1° 


Read  arrival  times 


Select  four  arrival  times  for  geometric  fix 


Calculate  four-station  fix 
{ 

i 

^  Calculate  velocities  from  velocity  function 


* 


Calculate  weights 

i 

Estimate  least -sous  ms  latitude  and  longitude 

i 

Compute  Lg  -  Lc  and  X  ^  -  \0 


differences  <0.1c 


- print  out  L0 ,  0 

4' 

Calculate  origin  times,  mean  origan  time, 
deviations,  standard  deviation,  and  Lapiacian 
of  variance. 


Calculate  Standard 
strength 


print  out 


power  levels  and  T-phase 


^  print  out 


Halt  or  proceed  to  next  event- 


The  computer  output  is  obtained  or  punched  cards,  printed  pates  arx’ 
on  a  plotted  map*  The  Printer  output  contains  details  of  the  solution 
process  -as  aids  in  reviewing  and  screening  the  results.  The  cards  for 
events  whose  solutions  appear  inadequate  are  pulled  from  the  deck*  If 
it  appears  that  the  results  can  he  improved  sufficiently  by  deleting 
inconsistent  input  data,  the  solution  is  recomputed. 

The  data  from  remaining  output  cards  is  then  listed  on  multilith 
masters.  A  sample  listing  is  shown  in  Figure  8.  The  explanation  of  the 
columns  is ; 

Columns  1-5  M  D  H  M  S  Greenwich  Mean  Time  of  T-phase  source 


Columns  6-7  LAT  LONG  Geodetic  coordinator  of  T-  phrse  source 


Column  8 


Name  descriptive  of  10°  square  contain¬ 
ing  T-phase  source 


Column  9 


Standard  deviation  of  origin  times 


Column  10 


CGNV 


Convergence  of  the  gradient  (Laplacian) 
of  variance  of  origin  times  in  seconds 
squared  per  malli radian  squared 


Column  11 


Num  r  of  hydrophones  read 


Columr  12 


T-phase  strength  in  decibels  relative 
to  0.1  microbar 


Lists  will  be  mailed  at  regular  intervals  to  interested  addressees. 
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The  Fortran  IV  compute j 
written  by  James  Sasser* 


program  comprised  of  7u0  punched 
The  mathematics  were  checked  by 


Roger  Norris*  William  M„  Adams  and  John 
of  Geophysics  and  David  Potter  of  General 


Northrop  of  Hawaii  Institute 
Motors  Defense  Research 


Laboratory  also  contributed  through  helpful  discussions..  This 


research  was  funded  by  the  Advanced  Research  Projects  Agency  through 
contract  Nonr  37ii6(01)  with  the  Office  of  Naval  Research. 
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